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Abstract: We study the possibility of emission of the saxion, a superpartner of the axion, 
from SN1987A. The fact that the observed neutrino pulse from SN1987A is in excellent 
agreement with the current theory of supernovae places a strong bound on the energy loss 
fS| ■ into any non-standard model channel, therefore enabling bounds to be placed on the decay 

constant, f a , of a light saxion. The low-energy coupling of the saxion, which couples at 
high energies to the QCD gauge field strength, is expected to be enhanced from QCD 
scaling, making it interesting to investigate if the saxion could place stronger bounds on 
^"'j f a than the axion itself. Moreover, since the properties of the saxion are determined by f a , 

a constraint on this parameter can be translated into a constraint on the supersymmetry 
breaking scale. We find that the bound on f a from saxion emission is comparable with 
^ . the one derived from axion emission due to a cancellation of leading-order terms in the 

soft-radiation expansion. 
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1. Introduction 

Searches for the neutron electric dipole moment (EDM) d n place strong constraints on the 
QCD vacuum angle #qcd [!]• The latest incarnation of these experiments gives an upper 
bound of \d n \ < 6.3 x 10~ 26 e cm which can be related [2] to the CP- violating parameter in 
QCD, 9 Q cd, to give 6 QCD < (1 - 10) x 10~ 10 . Bounds on 9 QCD coming from 199 Hg EDM 
experiments [3] are somewhat stronger, Oqcd < 1-5 x 10 -10 . However, the conversion of 
the bound on c?( 199 Hg) to a constraint on Oqcd is model-dependent and subject to large 
uncertainties. 

A possible explanation [4-8] for the smallness of #qcd is the introduction of a new 
anomalous global C/(1)pq symmetry which is spontaneously broken at the Peccei-Quinn 
(PQ) scale f a . Associated with this spontaneously broken symmetry is a new particle, called 
the axion [9, 10]. The axion is the pseudo Goldstone boson of the broken PQ symmetry. It 
develops a mass m a ~ m^f-wl fa due to the axial anomaly. 

The present constraints on f a are 10 9 GeV < f a < 10 12 GeV. The lower bound on the 
allowed window comes from SN1987A and corresponds to an axion whose coupling is such 
that it can be abundantly produced in the nascent neutron star and yet interacts weakly 
enough to allow emission on a similar time scale to the neutrinos. This axion provides a 
new efficient cooling mechanism and would lead to a shortening of the neutrino pulse from 
SN1987A, in conflict with observations, if f a is too small [11]. The upper limit on the 
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allowed window comes from cosmological arguments. It corresponds to axions so weakly 
coupled that they never thermalized in the early universe. At high temperatures the axion 
field value is misaligned with the value it takes at low temperatures and as the universe 
cools the axion begins to "roll down" towards its minimum [12]. This leads to coherent 
oscillations of the axion field. The requirement that the energy density in these oscillations 
does not overdose the universe leads to the upper bound on /„. 

In this paper we are interested in models which are supersymmetric and use C/(1)pq 
to solve the strong CP problem. If supersymmetry (SUSY) is combined with [7(1)pq there 
are superpartners to the axion, which are the fermionic axino and the bosonic saxion. For 
a given axion model the saxion's properties are determined by SUSY. In particular, for 
gauge mediated models [13,14], the saxion mass m s is comparable to the gravitino mass 
m 3/2 ~ 10~ 2 eV — 1 GeV. The saxion can be light and could be emitted from the interior 
of a supernova where the temperature is of order 20 MeV. 

We find that the saxion couples in a fashion similar to the dilaton [15]: It couples 
to the QCD gauge field strength in the ultraviolet (UV) and this coupling is enhanced 
in the infrared (IR) due to QCD scaling. This enhancement leads to the possibility that 
astrophysical bounds on f a from the saxion could be stronger than those from the axion. 
In certain gauge mediated models of SUSY breaking (see [14], for example) the PQ scale, 
f a , is derived from the SUSY breaking parameters. Any constraints on f a can be turned 
into constraints on these SUSY breaking parameters. 

Supernovae (SN), in particular the recent SN1987A, have been used frequently to ob- 
tain bounds on new physics using the energy-loss argument. Our current theory of SN 
explains the shape and duration of SN1987A's neutrino pulse and is in agreement with 
the observed data. If there exists a new channel that competes with the neutrinos and 
transports a comparable amount of energy from the interior of the SN then the current 
description of SN1987A's neutrino signal is significantly altered. A criterion for the max- 
imum possible emissivity of a novel energy loss process has been given by Raffelt [11]. It 
states that the luminosity of the new particle(s) can not exceed the neutrino luminosity 
of L v ~ 10 52 ergs/s. For a proto-neutron star as developed during SN1987A this criterion 
translates to an emissivity bound of e < 10 19 ergs/g/s for the new channel, the saxion in 
our case. Saxion parameters that lead to an emissivity greater than this are ruled out. 

Several calculations of axion emissivity from the process NN — > NNa have constrained 
f a [16-20]. Until recently, the NN scattering amplitude was derived from a one-pion- 
exchange (OPE) potential. However, as shown in [19], it is possible to relate the amplitude 
for axion emission to the on-shell NN scattering amplitude which can be derived from the 
experimentally measured NN scattering data. In contrast to the OPE calculation, this 
method gives a model-independent result. 

Naively one might expect that such a model independent calculation could be repeated 
for the saxion emission channel AW — ► NNs. This, however, is not the case due to a 
cancellation of the leading order terms in the expansion around the soft radiation limit. It 
is exactly these terms that make a model-independent analysis possible. 

In the non-relativistic limit the saxion — being a scalar particle — couples to the nu- 
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cleon mass. Hence there is no dipole radiation. Axions, on the other hand, couple to 
spin, which is not conserved in nucleon interactions and therefore they can be emitted by 
bremsstrahlung. Diagrammatically, the leading order poles corresponding to emission from 
an external leg sum in the case of axions and cancel in the case of saxions. 

For this reason we will calculate the amplitude for saxion emission at tree level in the 
SU(3) chiral Lagrangian. Using this Lagrangian to model the nucleon-nucleon interaction 
has several shortcomings, as explained below, and one might doubt the validity of our 
procedure. However, given the rather large uncertainties in the astrophysical parameters 
entering the calculation, it gives a sufficient estimate for the emissivity. 

We find that saxion emission is comparable to axion emission suggesting that the bound 
on f a from saxions could be at least as significant as that from axions. 2 However, several 
simplifying approximations were made in order to carry out the calculation. Once a model- 
independent calculation becomes feasible it would be interesting to repeat this calculation 
without relying on these approximations. Including the effects of this previously ignored 
emission channel could tighten the constraints on the PQ parameter space, if the UV 
physics is supersymmetric. 

This paper is laid out as follows. In Section || we derive the properties of the saxion at 
energies above the SUSY breaking scale. In Section || we compute the low energy effective 
coupling of the saxion. In Section || we discuss the soft radiation limit, calculate the saxion 
emissivity and derive bounds on f a . We conclude in Section |5[ 



2. (S) axions 

The gauge field piece of the standard model Lagrangian for a theory with a non-zero 
^QCD-angle is 3 

C = -^TrF^F^-^TrF^F^, (2.1) 

where 

Oqcb = e + Arg DetM (2.2) 

and g is the SU(3) gauge coupling constant. Classically, the bare 6 angle can be removed 
from the QCD Lagrangian by a chiral transformation of the quark fields. However, in the 
absence of a massless quark, this transformation is no longer a symmetry of the quantum 
theory since chiral transformations are anomalous. Therefore, 9 gets shifted by chiral phase 
rotations of the quark fields whereas Bqcd as defined in Eq. ( |2.2| ) is invariant. 

Consistency with experiment requires #qcd < 1CP 9 . Explaining the smallness of #qcd 
is the strong CP problem. The PQ solution to the strong CP problem postulates the 
existence of a new dynamical field, the axion a, which allows #qcd to vanish dynamically, 

1 This is similar to photon radiation from the collision of two charged particles with the same mass. In 
this case the center of mass coincides with the center of charge, there is no dipole radiation and no emitted 
bremsstrahlung. 

2 Here and throughout the paper we assume a light saxion, i.e. the mass of the saxion is less than the 
typical temperature in the neutron star (T ~ 10 — 50 MeV) and can be neglected. 
3 In the usual fashion, — e fl „ Q 0-F Q!/3 /2 and the trace is over group indices. 



-3- 



#qcd - * Qqcd + a{%(i)/fa- This is achieved by introducing an anomalous global symmetry 
U(1)pq spontaneously broken at the scale f a . 

Realizing U (1)pq in a supersymmetric theory necessitates putting the axion in a chiral 
supermultiplet. This multiplet is filled out by introducing its fermionic partner, the axino, 
and making the axion a the real part of a complex scalar field. The saxion s is the imaginary 
part of this scalar field and its couplings are uniquely determined by SUSY from those of 
the axion. 

In a supersymmetric theory the Yang-Mills action is most conveniently written in 
superspace notation as 

S = — Imr f d A xd 2 6TrW a W a , (2.3) 
8vr J 

where 9 is a two-component Grassmann variable and W a is a vector superfield which 
packages the gauge vector field A^, its superpartner the gaugino, and an auxiliary field. 
The coupling parameter r is chosen asr = ^ + ^ so that upon expanding Eq. ( p.3| ) 
one recovers the terms in Eq. ( |2.1| ). Making ^qcd dynamical corresponds to r -> r + 
where $ is a chiral superfield containing the axion, saxion, and axino. Expanding $ in 
terms of its component fields gives a shift of the bosonic piece: r — ► r + 2^J~( a + The 
saxion s couples to F 2 and the axion a couples to FF, as before. 

The axion is the pseudo Goldstone boson of C/(1)pq; it would be massless were it not 
for the axial anomaly. Without SUSY the axion gets a vev below the PQ scale and develops 
a 'mexican hat' potential. The radial mode becomes heavy (on the order of the symmetry 
breaking scale) and the angular mode remains light, becoming the axion. With unbroken 
SUSY the superpotential is holomorphic and this radial direction remains flat, ignoring 
the effects of the anomaly. There is now an extra light mode corresponding to oscillations 
in the radial direction: the saxion. With unbroken SUSY there are two degenerate light 
modes whose mass is set by the effects of the axial anomaly. Once SUSY is broken this 
nearly flat direction is lifted and the saxion acquires a mass, independent of the axion. The 
specifics of the SUSY breaking mechanism determine the form of the potential for <3? and 
thus the masses of the saxion and axion. Gauge mediated models predict viable properties 
for the saxion: Its mass can be light, of 0(lO~ 2 eV — lGeV), and f a fits naturally in the 
allowed window [13,14]. 

3. Low Energy Saxion Dynamics 

As we run down in energies from above the PQ scale, f a , passing the SUSY breaking scale, 
A^SUSY, to below the QCD scale, Aqcd> we pass through various effective theories that 
describe the relevant degrees of freedom at that energy. Summarizing what was mentioned 
in the last section (E ~ \f(p, where q is a typical momentum transfer): 

• E > f a , Mstjsy: U(1)pq unbroken, SUSY unbroken. 

• f a > E > Mstjsy: U(1)pq broken, SUSY unbroken, degenerate saxion and axion 
with mass generated by the axial anomaly. 
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• -^SUSY > E > Aqcd: SUSY broken, superpartners become heavy, nondegenerate 
massive saxion and axion, quarks and gauge fields are relevant degrees of freedom. 

• Aqcd > E: Mesons, baryons, saxion, and axion are relevant degrees of freedom. 
The Lagrangian for the saxion-axion system at f a is, 

C = \{B r s? + \(8 f a? - Via..) - + JL,±) TrF* 



16vr 2 



1 ^"qcd + yj TtFF + («* Pli ~ m i9iH + jd^aqri^Xq^j , (3.1) 



where the quarks q% have mass rm with the index i running over quark flavors. Furthermore, 
V(a, s) is undetermined and X is a model dependent flavor mixing matrix. 

In order to compute the couplings of the saxion to matter at the low energies typical 
for neutron stars we need to first integrate out the heavy quarks. We can then match 
this theory onto a low-energy Lagrangian containing hadrons as the relevant degrees of 
freedom. The saxion's interactions can be written in terms of the renormalization group 
invariant quantities rrnqiqi and d^D^, where the latter is the divergence of the scale current. 
This method of evaluating the divergence of the scale current in order to calculate matrix 
elements of the F 2 operator at low momentum transfer has been developed in [21-24] 
and recently applied to systems involving the low energy coupling of the dilaton [15] and 
quarkonium [25]. 

At the scale f a , where the Lagrangian in Eq. ([O]) is relevant, d^D^ is given by 

a 6 

d„D» = Tjt = ^TrF 2 + - jmjmiqm. (3.2) 

9 i=l 

Here we have ignored the contributions from the sF 2 and aFF terms since they are sup- 
pressed by a factor of l// compared to the terms in Eq. (|3.2|) . Furthermore , (3 = dg/d log fj, 
is the SU(3) beta function for the coupling g given to leading order by (3 = —bog 3 /l6ir 2 . 
Above the SUSY breaking scale 6o is given by &o = 9 — Nf for the Minimal Supersymmetric 
Standard Model. Below the SUSY breaking scale it is given by the standard model result 
bo = 11 — 2Nj/3, with Nf being the number of flavors. Furthermore, 7 mi = drrii/dlog fi 



is the mass anomalous dimension. Using Eq. (3.2) we can rewrite the saxion-gauge field 
piece of Eq. (|3.1|) as, 



(3.3) 



We can now run down in scale from f a to Aqcd via -^susYi matching across the quark 
mass thresholds as we go. Working to leading order in a s = ff 2 /47r and neglecting 7 m , 
which is small compared to unity, we arrive at 



Cs = lforV* ( d ^ ~ ^ ^l^j ' ( 3 - 4 ) 
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where the enhancement factor k is given by 



4 */ 9 ~\ og J^ + nA. (3.5) 



a,{U) 3 °GeV 

This running gives sizable enhancements to the saxion coupling. For instance, for f a = 
10 9 GeV, a SUSY breaking scale of 10 4 GeV, and A QC d = 200 MeV this QCD scaling 
results in an enhancement of k ~ 26 over the naive coupling. 



4. Calculating the Emissivity 

As detailed in the introduction, the Raffelt criterion limits any mechanism that competes 
with the neutrinos in removing energy from the neutron star to a maximum emissivity 
of 10 19 ergs/g/s. The emission of light saxions from the core of the neutron star would 
provide such an energy-loss mechanism. 

Emission rates have been calculated for the emission of other light particles, such 
as neutrinos (NN -> NNuv) [19,26], axions (NN -> NNa) [16-19], scalars (NN -> 
NNs) [27] , or Kaluza-Klein (KK) gravitons (NN — > NNg) which come about by introduc- 
ing large gravity-only extra dimensions [28-32]. There have also been several calculations 
for the emissivity of these particles where it has been assumed that the relevant degrees of 
freedom in the core are quarks and gluons, rather than nucleons [33,34]. 

In this work we consider a regime with nucleons as the relevant degrees of freedom. In 
the case of neutrinos, axions, and KK-gravitons it has been possible [19,30,31] by using low 
energy theorems [35,36] to relate the emission rates for soft emission to the on-shell nucleon- 
nucleon scattering amplitudes which can be extracted from NN scattering phase-shift data. 
Such a calculation has the advantage that it is model-independent, which distinguishes it 
from the often used one-pion-exchange approximation where the nucleon-nucleon amplitude 
is due only to a single pion exchange. A model-independent calculation is possible for the 
following reason: If the energy u of the radiated particle is low compared to the momenta 
of the nuclei (low radiation limit) then the main contribution to the emissivity comes from 
bremsstrahlung radiation where the radiated particle couples to an external nucleon line. 
These processes are 0(l/u>), they exhibit al/w pole due to an intermediate nucleon being 
nearly on-shell. The pole amplitude dominates over diagrams where the radiated particle 
couples to an internal line [which are O(oj )]. It is this dominance that renders the coupling 
of the emitted particle to unknown strong interaction vertices subleading and enables the 
model-independent calculation. 

As explained in Section |l|, a similar dominance of bremsstrahlung is absent for the 
saxion process NN — ► NNs. In each of the four bremsstrahlung-diagrams there is an 
infrared pole but these cancel in the sum — diagrams where the saxion radiates off one of 
the internal lines are not suppressed and low energy theorems cannot be constructed for 
saxion emission. 

One way of performing a model-independent calculation would be to use effective 
field theory (EFT) methods which have seen some advancement in recent years (for a 
review see, for example, [37]). Moreover, it has been shown by Beane et al. [38] that 
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an expansion of nuclear forces about the chiral limit is formally consistent and seems to 
converge as suggested by numerical evidence. Beane et al. investigated a toy theory of 
nucleons interacting with a potential consisting of three Yukawa exchanges (m p = 770 MeV, 
m a = 500 MeV, and = 140 MeV) where the masses and couplings where chosen in 
order to reproduce the scattering length and effective range in the 1 S'o nucleon-nucleon 
channel. By treating the pion part of the potential perturbatively (the pion coupling 
oc g\m^) this model shows convergence up to nucleon momenta \p\ ~ 250 MeV. While 
this progress is very encouraging, these methods are not yet developed far enough to be 
used at momenta > 300 MeV as are relevant for supernovae. Therefore we will calculate the 
amplitude for saxion emission in a SU(3) chiral Lagrangian for baryons and the pseudo 
Goldstone bosons associated with chiral symmetry breaking. OPE models, based on a 
SU(2) chiral Lagrangian have been used in the past for the case of neutrino, axion, and 
KK-graviton emission. Although the OPE captures the relevant physics of the nucleon- 
nucleon interaction in the long distance range, it fails to do so at short distances where 
the hardcore part of the potential becomes important. Although one may naively think 
that these contributions are unimportant for nucleon-nucleon interaction at non-relativistic 
energies, this is not true. Fine-tuning of the underlying parameters renders the short- 
distance contributions comparable to the long-distance ones. 

In cases where it is possible (e.g., axions and neutrinos) to calculate emission rates 
using both OPE and model-independent techniques, the two methods differ by about a 
factor of three due to the crude nature of the OPE assertion. We expect our calculation 
to differ from reality by a similar amount. This is sufficient for our purpose and will not 
change our conclusions considerably. 

Before proceeding we will briefly discuss other saxion emission processes and argue 
why AW — > NNs is dominant. A possible emission channel is the Compton-like process 
P7 — ► Ps, where P and 7 refer to a proton and photon, respectively. The leading order 
1/ll> piece in this process cancels between the two diagrams with uncrossed and crossed 
external proton lines and the process contributes at O(oj°). In addition, this process is 
suppressed by the fine-structure constant 1/137 and the proton number fraction Yp in the 
neutron star which is assumed to be small. Another emission channel, Nit — > Ns, where 
7T is either a thermal pion or a pion condensate, starts contributing at 0(u ) by the same 
argument. Also, thermal pions, being bosons, obey a Bose-Einstein distribution and their 
number density is suppressed by exp(-m 7r /T) rs 10 -3 . . . 10~ 2 for T « 20 ... 30 MeV. The 
existence of a pion condensate is rather uncertain and as such we will ignore its possible 
effects. In short, since the saxion- nucleon coupling is a scalar coupling it will commute 
with every other vertex and hence the 1/uj pole will always vanish in the soft-radiation 
limit. 

Before detailing our calculation for NN — > NNs, let us point out several effects we 
have left out in this paper. Throughout the calculation we use in-vacuum amplitudes to 
calculate processes in a dense medium. Whether this approximation is appropriate has 
been investigated, for axion emission, in Ref. [39]. It turns out, that high-density effects 
are significant and can modify the emission rate by a factor of about an order of magnitude. 

Next, the emission rate can be suppressed significantly due to multiple-scattering ef- 
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fects as pointed out by Raffelt and Seckel [40]. This effect, which is analogous to the 
Landau-Pomeranchuk-Migdal (LPM) effect for electromagnetic bremsstrahlung by rela- 
tivistic electrons, applies if the nucleon collision rate exceeds the oscillation frequency of 
the emitted radiation. If the time to emit radiation is greater than the time between colli- 
sions in the medium then the nucleons will undergo multiple collisions before emission. It 
is the interference between these collisions that will suppress the radiation. One incorpo- 
rates this effect be assigning the nucleon a finite decay width, T. This can be done by the 
following substitution in the squared matrix element: 



(4.1) 



UJ 2 uj' x 



+ r 2 /4' 



A typical width is given by [41] 



T = aT 2 



1 + 



V2ttT/ 



(4.2) 



where a « 0.2 MeV" 1 . The LPM effect sets in when r ps lo, or T > 10 MeV. It suppresses 
the emission at high temperature, making the late-time cooling phase of the star more 
important. 

For radiation of particles that couple to the spin of the nucleon, such as axions or 
neutrinos, there is an additional suppression. If the correlation length of the nucleons 
is small compared to the formation length of the emitted radiation, then the emitted 
particle couples to an averaged spin since the nucleon-nucleon interaction can flip spin. 
This averaged spin goes to zero as the nucleon correlation length decreases. For axion 
emission this effect reduces the limit on f a by about a factor of 2 [42]. The saxion, however, 
is a scalar and does not couple to the nucleon spin. We therefore expect no suppression of 
the emissivity due to this effect. 

4.1 SU(3) Chiral Lagrangian and Saxion Coupling 

For low-momentum processes chiral-perturbation theory is the suitable framework to work 
in. In contrast to OPE calculations, which are based on SU(2), we use in our calcu- 
lation a SU(3) chiral Lagrangian. The nucleons in the core of the neutron star have a 
typical momentum of ~ 300 MeV which suggests the inclusion of the strange quark in 
the calculation. Moreover, as will be explained in more detail below, the strength of the 
saxion-nucleon coupling depends on the masses of the octet baryons, which necessitates 
SU(3). As a consequence, at tree-level there are not only 7r° but also r\ exchanges between 
nucleons. 

The strong dynamics of the octets of pseudo Goldstone bosons M and baryon octet B 
is described, to lowest order, by the Lagrangian 

f 2 

+Tr£ t (iD + —\ B + DTrB^G i {A i ,B} + FTrB^a i [A i ,B] (4.3) 
y m-B I 
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where = dp + Vp and vtiq is the quark mass matrix ffiq — diag(?ra u , m,d,m s ). Furthermore 



£ = exp^— j, £ = ^, M = -AV 

^ = \ + £c^) , A^= l - (id^ - gdpZ) , (4.4) 

and the \ a are the SU{3) generators. We use the convention where the decay constant 
/ = 132 MeV. 

The proton number fraction Y p is typically small in SN, especially during later times in 
the SN's evolution when most of the protons have neutronized. Therefore, and in order to 
unclutter our formulae, we specialize our treatment at this point to ran scattering by taking 
Y p = 0. Although we could carry out the calculation for the general case of a proto-neutron 
plasma specified by Y p , we expect the inclusions of protons to make little difference in our 
result. 

To lowest order in the chiral expansion, the contributing diagrams are those where a 
single 7T° or r\ is being exchanged between two neutron lines. Therefore, we need to just 
keep neutrons n as well as 7r° and rj. Then the Lagrangian in Eq. Q4.3p reduces to 

V 2 

C = n^(ido H — 7i~) n ~ dn^a ■ Vtt° n — g v n^a ■ V77 n 
£iv± n 

+~ VW - ^(vr ) 2 + IdprPv - ^V, (4.5) 

where n is the neutron 2-spinor and M n is the mass of the neutron. The couplings of 
neutrons to 7r° and rj are 

D + F D-3F 

Note that the Vp coupling has disappeared as it involves a baryon-baryon-pion-pion cou- 
pling which does not contribute to saxion emission at tree-level. The meson masses can be 
calculated in terms of the quark masses, 



4A, > 2 4Xm u + m d + im s 
m n = j^{m u + m d ), m v = . (I.<) 



The parameters D and F can be determined from fitting to the measured hyperon 
semileptonic decays. This has been done by several authors [43-47]. At tree level a typical 
fit gives D = 0.79 ± 0.10 and F = 0.47 ± 0.07, where the errors are highly correlated, 
resulting in the combinations D + F = 1.26 ± 0.08 and 3F - D = 0.65 ± 0.21 [45]. If 
one includes chiral loop corrections, which are 0(m s lnm s ), a typical best fit gives D = 
0.64 db 0.06 and F = 0.34 db 0.04 [45], where intermediate decuplet baryons have been 
included. For our tree level calculation we use typical values of D = 0.8 and F = 0.5. 

To calculate the coupling of saxions to n, ir° and rj from the Lagrangian in Eq. ( fTjj ) 
we first need to calculate the trace of the energy-momentum tensor Tp V from Eq. ([4.3;). 
Making use of the equation of motion for neutrons, 4 we find 

djj^ = T£ = M n nn + 2m 2 ^° 2 - (^ t vr ) 2 + 2m 2 r ? 2 - (d^f . (4.8) 

4 This is justified since we only calculate diagrams where the nucleons are external states and therefore 
on-shell or nearly on-shell. 
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Since the saxion couples to d^D^ — Ya=i m iQiQi we need to calculate the matrix elements 
of this operator sandwiched between our low-energy degrees of freedom, neutrons, ir°, and 
rj. Matching 

3 22 
£ m,q m — -ATr(m 9 St + m,S) = -^(vr ) 2 + -frf (4.9) 

i=l 



we get 



and 



{ir°{q')\d^ - mqiqi^iq)) = ~2q' • q + 3m 2 (4.10) 

8=1 



MgOlfyfl" " E ™iqi<H\v(Q)) = -V • g + 3m 2 ,. (4.11) 

i=l 

Similarly we obtain for the neutrons n 

3 

<n(p)|5 M ^ -X)migigi|n(g)) = M (4.12) 

2=1 

where M = 687 ± 75 MeV [15, 48] and can be calculated from the masses of the octet 
baryons and the pion-nucleon S-term. The error for M stems from an estimated 30% error 
typical for SU (3) violation. Finally, the couplings of the saxion to the nucleon-pion and 
nucleon-eta vertices vanish, 

3 3 

i=i i=l 
Putting these interactions into the saxion Lagrangian [Eq. Q3.4j )] we arrive at 

Cs = lbW [Msrfn + ^m 2 S (7r ) 2 - s(d^°f + ^m 2 srj 2 - s(d^) 2 ^j . (4.14) 

Before calculating the saxion emission amplitude in detail using the chiral Lagrangian 
and saxion coupling given above, we will do a comparison between the cross sections for 
saxions and axions in nuclear matter based on dimensional arguments. 

4.2 Comparison between Saxion and Axion Emissivity 

The coupling of axions to nucleons is given by 

—N^a-VaN (4.15) 

J a 

so that the lowest order in the uj expansion of the cross section for NN — > NNa goes like 

1 uj 2 (M n \ 2 
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k k 
k k 

Pi ^ . fc Ps 

Pi ». T » ■ » P3 Pl *- ? «■ P3 T 

i i 1 

i i i — ► — ^ 

i i i 

i i i 

P2 » • * £>4 P2 »* • P 4 1 



£>2 P4 

(a) (b) (c) 

Figure 1: The three types of diagrams for the process nn — > nns. Including crossings, there 
are 16 diagrams (8 pion and 8 eta exchanges) of the type shown in panel (a), 4 diagrams of type 
(c), and no contribution from type (b) diagrams. A solid circle denotes a meson- nucleon-nucleon 
vertex, while a solid (empty) square denotes the coupling of the saxion to a nucleon (meson) line. 
The saxion coupling to the meson-nucleon-nucleon vertex is denoted by an empty circle. A dashed 
line denotes either n° or 77. 

Here, the (M n /u) 2 term is due to the l/u> pole coming from a nearly on-shell intermediate 
nucleon. Using the Lagrangian for the coupling of the saxion to nucleons given in Eq. (|3.3| ) 
we can write a similar cross section for NN — > NNs: 

" ~ P?M*5^- <4 ' 17) 

Now, because of the aforementioned 1/uj pole cancellation in the case of the saxion, the 
piece of lowest order in u> is 0(u> ). Comparing the expected saxion and axion emissivities 
we find 

e a n 2 M 2 

1(T 2 . (4.18) 



e a 256vr 4 M2 

So on dimensional grounds we expect the saxion emissivity to be suppressed be about 
2 orders of magnitude compared to the axion emissivity. 

4.3 Matrix Element 

The relevant saxion emission processes are shown in Fig. [l]. There is a total of 16 bremsstrahlung 
type diagrams (4 direct and 4 crossed pion exchanges, 8 eta exchanges) [referred to as type 
(a) diagrams]. Moreover, there are 4 diagrams where the saxion couples to the exchanged 
meson [type (c)]. The type (b) diagrams do not contribute. 

As an example, the diagram pictured in Fig. |l](a), where the exchanged meson is the 
7T°, gives a contribution of 



16vr2/a 



(P4-P2)i{P4-P2)j 



1 {M + E 3 

- + 1 



(P4 -P2) 2 ~ ml 



n 



i a i nxn\a j n 2 . (4.19) 



2E Z \ uj 

Similarly, the diagram depicted in Fig. |l|(c) for the 7r° exchange contributes 
~ iK 9l r \ r \ 1 1 

, R 27 IPl -P3)i{P4 -Pl)r< \2 2 7 \2 2 

16n 2 f a (Pi - P3) - (p 2 - P4) 2 - 

x [3ml - 2(pi - p 3 )(p4 - p 2 )) nl^nx n\a J n 2 {4. 20) 
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Here, o"j are the Pauli matrices and E{ = M + p~i 2 /2M in the non-relativistic approx- 



imation. Note that in Eq. (4.19) we need to keep terms up to O(lu ). These will become 



the main contribution from the type (a) diagrams as the pole drops out. 

The squared amplitude can then be calculated after summing over spins in terms of 
the 4-momenta k = p 2 — Pa and I = P2 — P3 and is found to be 



\ M \ 2 = T,\M ia) +M ib) +M (c) 

spins 
K 2 M 4 



{(K x + K 2 ) 2 k 4 + (Lx + L 2 ) 2 l A - (Ki + K 2 ){L X + L 2 ){2[k ■ I] 2 - k 2 l^l) 



where 



A + (4.22) 



Eqm \k 2 — m 2 k 2 — m 2 

^ = 9 2 n (3ml-2k 2 ) + 9 2 (3m 2 -2k 2 ) 
(k 2 — m 2 ) 2 (k 2 — m 2 ) 2 

and L\ and L 2 are defined analogously with k — > I. 

In this calculation we have assumed the mass of the saxion m s to be small compared 
to the nucleon momenta and have neglected it. 

4.4 Phase-space integration 

In order to calculate the emissivity due to saxion s radiation from the two-body scattering 
reaction a+b — > a+b+s in neutron matter, where the labels a and b represent a neutron, we 
have to integrate over the phase space of incoming (pi, p 2 ) and outgoing (ph, pa) nucleons 
as well as that of the emitted saxion (k s ). The formula for the emissivity is 



d£x 
dt 



x(2it) 4 5Hpi+P2-P3-P4 ~ k)Suo\M\ 2 (4.24) 



where the functions fa are Pauli blocking factors given by /j = (exp[(Ei — p)/T] + 1) _1 and 
S = 1/4 is the symmetry factor for nn scattering. 

In the soft-radiation limit we can neglect the saxion momentum k s in the momentum- 
conserving delta function. Furthermore, by introducing the total momentum P = p~[ +p 2 = 
p~3 + p~a we can reduce the number of integrations by exploiting spherical symmetry and 
momentum conservation to obtain 

dS 1 /*oo /'OO /'OO ^1 rl pI'K 

— — = -T^: / dP\P\ 2 / dk\k\ 2 / dl\l\ 2 / dcos9 k / dcosOi / d<f> 
dt 2 ib 7T 8 J Q J J J_ x J_ x J Q 



x f hf 2 (l - / 3 )(1 - U)\M\ 2 . (4.25) 

The angles Ok and 9\ are defined by k ■ P = \k\\P\ cos 9^ and I ■ P = \l\\P\ cos 6i, while <j) is 
the angle between the projections of k and I onto the P-plane and is given by 

cos 6*cm = cos <p sin 9^ sin 9i + cos 9^ cos 9i (4.26) 



- 12 - 



where cos #cm is the center of momentum (COM) scattering angle. Also, uj is constrained 
by the energy delta function to be uj = E\ + E 2 — E3 — E4. 

The expression in Eq. ( 4.25| ) is valid for all temperatures and chemical potentials. In 
the general regime it must be evaluated numerically. 

The result simplifies significantly if one assumes either a degenerate or a highly non- 
degenerate neutron gas in the star. Although nuclear matter at densities of a few times 
nuclear matter density and a temperature of 10 — 50 MeV is neither degenerate nor non- 
degenerate, it is nevertheless instructive to calculate saxion emissivity in these two limits. 

4.5 Limiting Case: Degenerate Neutron Gas 

For a degenerate neutron gas one can assume that saxion radiation only arises from scat- 
tering involving neutrons near the Fermi surface in the initial and final states. Assuming 
soft radiation the momentum \pl\ of the neutrons in Eq. ( [4.24j ) can be set to pf causing a 
decoupling of the energy and angular integrations: 

= 2^10 r duJUj2 / dE i dE 2dE 3 dE 4 5(E 1 + E 2 -E 3 -E 4 - u/)/i/ 2 (l - / 3 )(1 - h) 

x J dnxdn 2 dn 3 dn 4 5 3 (pi + f 2 - f 3 - pa)\m\ 2 (4.27) 

where the 5Qi is the angular integration of the pi. 

Carrying out the energy [49,50] and angular integrations one finally integrates over co 
to find 

H = ¥^ PFT6 [ da ^J^ecmlMf, (4.28) 

where cos a = P ■ p~\/\P\\pi\. The COM scattering angle 6 C m is the angle between the 
projections of pi and p~s into the P plane. 

4.6 Limiting Case: Non-Degenerate Neutron Gas 

In the limit where the system can be treated as a non-degenerate system the initial-state 
distributions become (non-relativistic) Maxwell-Boltzmann distributions, 

/,(«) = f(^) 3 %(-H) *- « = <«*>) 

where 

n=2 Sw? m) (43o) 

is the neutron number density and the factor of 2 is for the two spin states of neutrons. 
The final-state blocking factors can be set to unity since basically the complete final-state 
phase space is available. Then, starting from Eq. ( [4.25 ), the integration over P can be 
done analytically, leaving 

dtX = n 3/2 15/2 

dt 2 14 7T 9 / 2 

x |fc| 4 |Zj 4 cos 2 ^exp I - ' ri ' j:"" 1 """"" I \M\\ (4.31) 
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Figure 2: Saxion emissivity for fixed f a = 1 GcV, k = 1 at nuclear matter density depending 
on temperature. The crosses are values from the full phase-space integration, the solid and dashed 
lines denote the result from the degenerate and non-degenerate matter calculations. 



4.7 Numerical Results 

We show in Fig. ^ the emissivity for nuclear density n = no, where uq = 0.17fm~ 3 = 
1.36 • 10 6 MeV 3 , and varying temperature including the limiting cases of degenerate and 
non-degenerate matter. In the figure we have set f a = 1 GeV and k = 1. As expected, 
the degenerate approximation is good at low temperature while the non-degenerate one is 
better suited to a high temperature regime. 

For a typical temperature of T = 20 MeV and n = uq the emissivity is given by 



^f^Vs.SxlO 32 ^ (4.32) 
V fa J g-S 



so that, upon using Raffelt's criterion of 10 19 ergs/g/s and Eq. ( |3.5| ), we can solve for the 
PQ scale f a which gives a bound of 

/„ > 1.5 x 10 8 GeV (4.33) 

that is close to the bound of 10 9 GeV coming from the axion. However, if the temperature 
in the neutron star is as high as 50 MeV the same calculation gives a bound as strong as 
that from axion emission. 
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The inclusion of the r], typically absent in OPE calculations, increases the emissivity 
by only about 10% because it is four times as heavy as the pion. Furthermore, we expect 
the LPM effect to weaken our bound for T = 20 MeV by a factor of ~ 5. 

5. Conclusions 

If both PQ symmetry and SUSY are realized in nature, a new scalar particle, called the 
saxion, is introduced as a supersymmetric partner of the axion. The properties of the 
saxion are determined uniquely by SUSY and linked to those of the axion. Both particles, 
if they exist and are light, are produced in SN. In the past, several calculations using the 
energy-loss argument for the axion have been made in order to place bounds on the PQ 
scale f a - In this paper we used saxion emission to bound /„. 

The fact that the axion couples to TtFF while the saxion couples to TrF 2 at an 
energy scale of f a causes the two particles to couple very differently to nucleons at energies 
typical to neutron stars. In the limit that there is a massless quark the strong CP problem 
disappears and there is no need for an axion. Hence, for a quark mass approaching zero 
the axion-nucleon coupling goes to zero. In the same limit, however, the saxion does not 
decouple since it couples to the QCD field strength. One might therefore expect that the 
saxion coupling is enhanced over the axion coupling. There is also the added effect of 
enhancement due to QCD scaling in the case of the saxion. On the other hand, saxion 
emission is suppressed because the 1/lo pole, which is present for the axion, is absent. Thus 
it is not clear a priori how the bounds on f a coming from the two particles compare. 

Ideally one would like to repeat the model-independent calculation of the axion emis- 
sivity [19] for the saxion. This turns out not to be possible because the 1/lo pole is not 
present. Therefore, low energy theorems, which enable one to relate the emissivity to the 
on-shell nucleon-nucleon scattering amplitudes are no longer applicable and one is forced 
to confront the full details of the nucleon-nucleon strong interaction. 

In the past the nucleon-nucleon potential has frequently been approximated by an 
OPE. This OPE approximation is of a crude nature since it misses short-distance physics 
which is more important than one might naively expect. Therefore its validity is difficult 
to quantify. However, in the case of neutrino and axion emission comparisons of the OPE 
approximation with a model- independent calculation have shown that, surprisingly, the 
two methods give similar results to within about a factor of three. In our calculation of 
saxion emission we used an SU(3) chiral Lagrangian, calculating the cross section to tree- 
level. Since this method is closely related to the OPE approximation, one might hope that 
it does equally well. 

Combining our result for saxion emission with the Raffelt criterion we find a bound 
on f a that is very close to the one coming from axion emission. In view of the significant 
amount of uncertainty in our calculation this raises the possibility that saxion emission 
could significantly raise the lower bound on f a . We hope that the exciting prospect of 
further tightening the bound on f a gives motivation to a model-independent repetition of 
our calculation once the necessary EFT tools become available. 
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